Casimir force on interacting Bose-Einstein condensate 
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We have presented an analytic theory for the Casimir force on a Bose-Einstein condensate (BEC) 
which is confined between two parallel plates. We have considered Dirichlet boundary conditions for 
the condensate wave function as well as for the phonon field. We have shown that, the condensate 
wave function (which obeys the Gross-Pitaevskii equation) is responsible for the mean field part of 
Casimir force, which usually dominates over the quantum (fluctuations) part of the Casimir force. 

PACS numbers: 03.75.Hh, 03.75.-b, 42.50.Lc, 05.30.Jp 



I. INTRODUCTION 

Confinement of the vacuum fluctuations of the elec- 
tromagnetic field between two plates give rise to a long 
ranged Casimir force [IJ. The experimental verification 
of the Casimir effect is a confirmatory test of quan- 
tum field theory. Although the theory of Casimir force 
was first given [l[ for zero temperature (T), yet it can 
be generalized for any range of temperature (a, @ and, 
for any dielectric substance 0, [1[ between two dielec- 
tric plates. In general, Casimir forces are always present 
in nature when a medium with long ranged fluctuations 
is confined to restricted geometries Consequently, 
it has been generalized for thermodynamical and criti- 
cal systems jol— Tl~2l] . Recent measurements [l3l - [l6| of the 
Casimir forces for these kind of systems have drawn in- 
terests to the theoreticians 

@, El EMI- Although the 
Casimir force for a realistic BEC has not been measured, 
yet the Casimir-Polder force for the same was measured 
by the experimentalists of the Ref . [2(} ■ A successful the- 
ory for the Casimir-Polder force was also obtained by 
the authors of the Rcf.]21]. Since many new interesting 
phenomena are being discovered with the experimental 
studies [22-24] of BEC, we are optimistic to have the 
measurement of the Casimir force on the BEC. On this 
issue, the Casimir effect on BEC has been the subject of 
a number of theoretical works within the last few years 




One of us recently studied [H, [2!| the temperature 
dependence of the Casimir force for an ideal Bose gas 
below and above of its condensation temperature (T c ). 
Below the condensation point, it vanishes as a power 
law of the temperature and is long ranged. High above 
the condensation point, it vanishes exponentially and is 
short ranged. The role of geometry in the calculation of 
Casimir force for trapped ideal Bose gas was also explored 
in the Ref. c 29] . Different geometry causes different power 
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law behavior [H [29[. The Casimir force for the non- 
interacting BEC was also known to be vanished at the 
zero temperature 0, [2l|. But, the Casimir force for in- 
teracting BEC at zero temperature may not vanish. The 
same force due to zero-temperature quantum fluctuations 
(phononic excitations) in an interacting one dimensional 
harmonically trapped (but homogeneous!) BEC was also 
obtained recently [30| . It is to be mentioned, that, the 
theory of the Casimir force for a realistic (3D harmoni- 
cally trapped and interacting) BEC has not been given 
so far. 

We already have mentioned, that, the Casimir force 
depends on the geometry of the system. From this de- 
pendence, we can expect, that, the different geometry 
of the same system may not change the basic nature of 
the Casimir force. Hence, instead of the studying the 
Casimir force for a realistic BEC, it is relevant to study 
the Casimir force for a 3D interacting BEC in the plate 
geometry. On this issue, the Casimir effect due to the 
zero-temperature quantum fluctuations (phononic exci- 
tations) of a homogeneous weakly-interacting dilute BEC 
confined to a parallel plate geometry, has recently been 
studied with periodic boundary conditions [13, H2[ . 

Within the plate geometry, the BEC is considered to 
be homogeneous when the plate separation is infinitely 
large. Otherwise, the condensate is inhomogeneous, and 
it obeys Dirichlet boundary conditions. Since a macro- 
scopic number of particles below T c occupy the single 
particle ground state, the theoretical predictions of the 
BEC for T — » 0, match well with the experimental results 
which are actually obtained for T S, T c . Hence, we are 
interested in calculating the Casimir force for T on 
a 3D interacting inhomogeneous BEC confined between 
two parallel plates. 

Our calculation starts from the standard grand canon- 
ical Hamiltonian (H) for the interacting Bose particles. 
Then we obtain the grand potential for the condensate, 
and quadratic (quantum) fluctuations for the phonon 
field. The grand potential is the mean field part of 
the H. The condensate wave function obeys the Gross- 
Pitaevskii (G-P) equation [33. From the consideration 
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of the Dirichlet boundary conditions, the condensate be- 
comes inhomogeneous. Then we obtain the mean field 
force acting on the plates in terms of the maximum of 
the inhomogeneous condensate wave function. The max- 
imum of the condensate wave function is asymptotically 
obtained for very small and very large limits of the plate 
separation. For the whole range of the plate separation, 
we obtain an interpolation formula, which matches very 
well with the asymptotic solutions. We also justify this 
interpolation from the exact graphical solutions. We ob- 
tain the mean field Casimir force from the interpolation 
formula. Then we obtain the Casimir force for the vac- 
uum fluctuations of the phonon field with the Dirich- 
let boundary conditions. Finally, we compare the mean 
field and the quantum fluctuations contributions to the 
Casimir force. 



II. ELEMENTARY EXCITATIONS OVER A BEC 
IN THE PLATE GEOMETRY 

Let us consider a Bose gas of N identical particles to 
be confined between two infinitely large parallel plates of 
area A along the x — y plane and their separation along 
the z direction be L. For T — > 0, almost all the particles 
(Nq fa N) form the condensate, and the elementary ex- 
citations over the condensate be perturbatively treated. 
The elementary excitations over an unbounded conden- 
sate are obtained by the standard textbook approach 
[33l [HJ . We follow the same approach for obtaining the 
elementary excitations over a bounded condensate in the 
plate geometry. 

Let the position vector and mass of a single Bose par- 
ticle be denoted as r = r ± + zk and m respectively. The 
grand canonical Hamiltonian operator for such a system 
of interacting Bose gas is given by [33[ 



U = j # f (r) 



2m 



V 2 -/x *(r)d d r 



+ - J J I ^(r)^(r')V(T-r')^(T)^(r')\d 3 rd 3 r', 

(1) 

where fi is the chemical potential, V(r — r') is the inter- 
particle interaction potential, and ^(r) is the field op- 
erator for the Bose particles. For the simplest case, the 
interaction potential can be considered as V(r — r') = 
<7<5 3 (r — r'), where g is the coupling constant. The con- 
nection of this coupling constant with the s-wave scatter- 



ing length (o s ) is [33J g 



inha s 



The field operator in 



terms of the single particle orthonormal wave functions 
{(f>i(r)} is expressed as \&(r) = ^°^ o 0i(r)aj, where en 

and a] annihilates and creates respectively a Bose parti- 
cle at the state (r) . Within the perturbative approach, 
the grand canonical Hamiltonian in terms of the excita- 
tions S^>(r) (= ^(r) — V No(po(r)) be effectively written 



up to the quadratic order as [34 



n = Qo + Js¥(t)(- |^V 2 ^(r)d 3 r 

+ (J*(r)<y*(r))d 3 r, (2) 



where 

n = N 



2m 



|V</» (r)| 2 - M |^ (r)| 2 



gN 



l<Mr)| 4 K r 



(3) 



is the grand potential for the condensate, and where 
7Vo0o(r) terms and gn in the quadratic fluctuations are 
approximated by the bulk density (n) and the chem- 
ical potential (/x) respectively [34|. Within the per- 
turbative approach we can write the quantum fluctua- 
tions for the Dirichlet boundary conditions as S^(r) = 



where a T 



annihilates a Bose particle of x — y momentum pj^ = 

pj + Py j and energy & + sgf = £. The % 
in the Eqn.(2) can now be diagonalized in terms of 
the phononic operators through the Bogoliubov trans- 



formations I35 |j a , 
7T + 



= u 



■j ''p 



v 



pu 

v±,3"p±,j ' "pij"-p 

f p 2 /2m+gn n 1 / 2 _ ( p 2 /2m+gn 

V 2f.fa, .i) 2) ' U P±J 

Si 



P±,3 b -p ± ,j and 



j , where [33jj 



2e(px,i) 



IN 1/2 
2J ' 



and 



e(p±,j) = ( — [P 2 



9 n (j2 , 7r 2 /t 2 j 2 
1? 



„2 , K 2 h 2 j 2 N N 1/2 



Amgn 



(4) 



With the above transformations, we recast the Eqn.(2) 
in terms of the phononic excitations as 



1 V"^ ( , ,p\ 7T 2 ?l 2 j 2 N 

•2m 2mL 2 ' 



gn 



e (P-L^)&Pi,A_L,r 



(5) 



For T — > 0, there would be no phonon, and consequently, 
the grand canonical energy (< % >) of the system for the 
vacuum of the phonon can be obtained from the Eqn.(5) 
as 



£ = il + s (L,n), 



(6) 



where e (L,n) =iE Pilj j) " (f£ + l2rx) ~ 9 n ) 

is the contribution to the grand potential {£q) due to 
the quantum (vacuum) fluctuations of the phonon field. 
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Roberts and Pomeau [35] also correctly predicted this 
term from the original work of Lee, Huang and Yang 36] 
on a homogeneous condensate. 

From the second term of the Eqn.(6), we will obtain 
an expression for the Casimir force due to the vacuum 
fluctuations. Irrespective of the form of the second term, 
our expectation for the Casimir force due to the vac- 
uum fluctuations of the phonon field, would be similar to 
that due to the vacuum fluctuations of photon (electro- 
magnetic) field. Since, the Casimir force for the vacuum 
fluctuations of photon field is C ptn — — 2aol^ Q 



our 



expectation for that due to vacuum fluctuations of the 
phonon field (6&) would primarily be [III 0j5[ 



C, 



phn 



Air 2 hv(n) 
480L 4 ' 



(7) 



where the speed (c) of light is replaced by the speed 
(v(n) = \fgnjm) of phonon, and a factor 1/2 appears 
for a single polarization of a phonon in a BEC. 

Before going into the details of the quantum fluctu- 
ations part, let us calculate the mean field part of the 
Casimir force from the first term of the Eqn.(6). 



III. MEAN FIELD FORCE ACTING ON THE 
PLATES 

In the first term (fio) of the Eqn.(6), we can replace 
</>o(r) by 4>o(z) due to the fact, that, the condensate is 
placed between two infinitely large parallel plates, and 
the motion of the condensate is relevant only in the z 
direction. This term governs the equation of motion 



2m dz 2 



At + gN \(j)o(z)\ 2 )Mz 



(8) 



which is known as time independent Gross-Pitaevskii (G- 
P) equation. In the above analysis, (j)o( z ) has been con- 
sidered to be a complex function which can obviously be 
written in terms of the condensate density (tiq(z)) and a 

real phase (9(z)) as o (z) = ^&e i9 W. 

We are considering the fact that no macroscopic part 
of the condensate is moving as a whole within the plate 
geometry. This consideration allows us to write J — > 0, 
where J = n(z)^V9(z) is the particle current density. If 
J — > 0, we can ignore the phase part of the kinetic en- 
ergy density (^[(V v / no(z)) 2 +n (z)(V6l(z)) 2 ]), and can 



write 



i_)(vv^R) 



which of course allows the 



no current case to be applied. Under the same condition, 
we consider <fio(z) to be a real function for the rest of this 
paper. 

Now, with the substitutions a — 2/im/fi 2 and b = 



2mgNo/h 2 , the Eqns.(3) and (8) can be recast as 



and 



d 2 Mz) 
dz 2 



a<f> (z) + b<pl(z)=0 



(9) 



(10) 



respectively. 

From the Eqn.(9), we get the mean field force act- 
ing on the plates as F m f(L) = —^-\ : 



NAh 



( 1 ( d<Po \ 
\2 \ dz > 



|^(L) + U*(L)), where f£ 



dOo I 
dL \a 

on 1 _ 



t>l$fe)dz is associated 



NAh 2 ( L (d4o jl_td4a_ \ _ nr u dfo 

m JO y dz dL \ dz I Ul ^ dL 1 "^O dL 

with the change in the wave function with respect to 
the change in the plate separation, and the other part 
of ^§ L a -\ a is associated with the change in the integra- 
tion limit. Now, from the Dirichlet boundary conditions 
(</>o(0) = <f>o(L) — 0) and from the Eqn.(10), we recast 
the above expression of the mean field force as 



F (L)= NAh2 ( d ^ d ^ 
m \ dz dL 



2 I dz 



z=L 



.(11) 



Since 4>o(z) obeys Dirichlet boundary conditions, we 
can write it in terms of the Fourier modes as (j)o{ z ) — 
Y^jLifi(j) x ^{^jr)^ where fi(j) is to be determined 
from the Eqn.(10). However, we will calculate the 
mean field force without determining fi(j). From the 



Fourier expansion, it is easy to check that 



dL 







-^7 With these properties, the 



Eqn.(ll) becomes [1 



Fmf(L) 



NAh 2 1 i uip 
dz 



(12) 



This is to be mentioned that, the Eqn.(12) can also be 
obtained from the stress tensor [37[ of the G-P equa- 
tion. In the Ref . [37i] . the stress tensor was calculated for 
non-interacting bosons and fermions. The stress tensor 
(whose average over the statistical distribution is actu- 
ally the pressure at any point of the system) involves all 
possible orthonormal single particle states and their oc- 
cupations. If the inter-particle interaction be introduced 
in the formalism of the Ref. (37| , the Schrodinger equation 
for the ground state will become the G-P equation. With 
these considerations, the stress tensor of the Ref. [13] will 
give exactly the same force at the boundary (z = L) of 
the system as we have obtained in the Eqn.(12). 

Besides the Dirichlet boundary conditions </>o(0) = 
4>o{L) = 0, <f>o(z) is obviously symmetric about z = 4, 
and, since 4>o(z) is the ground state of the system, it must 



have a single maximum at z 



Let us now multiply 



the Eqn.(10) by 2^£f and integrate with respect to z from 
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z = j% to z = z. From this integration and with the fact 
that r = 0, we have 

(13) 

where $j(|f) * s an integrating constant. From the 
Eqn.(13) and from the Dirichlet boundary conditions, we 
recast the Eqn.(12) as [l9j 



Fmf(L) 



NAh 2 lf l9 fL 

—2 rn 2 

iVAfr 2 a 2 , 

?7(1 - r?), 



6 , 4 /£ 



6m 



(14) 



where rj — j-«^o('§)- From the Eqn.(14) we see that, we 
need only to know 77 as well as the maximum of (f>o(z) for 
obtaining the mean field force acting on the plates. In 
the next section we will write a transcendental equation 
for the maximum of (j>o(z). 



IV. ANALYSIS OF THE GROSS-PITAEVSKII 
EQUATION 

Eqn.(13) can once again be written in the differential 
form 



dz 



#o(z) 



^a{^)-4>l(z)) + \{^)-^)) 



(15) 



Now, integrating the Eqn.(15) from z — to z — z, we 
get 



1 



(16) 



v/«(l - ^ ^ v/(l-t 2 )(l-Mi 2 ) 

^ jai- From the Eqn.(16), 



where t = and M = - JL 

we get the ground state [3 

f2M 



2EllipticK[M] 



x JacobiSN[2EllipticK [M]-,M], (17) 

where EllipticK[M] = § [l+f- + ^M 2 +^M 3 +...] is the 
complete elliptic integral of first kind, and JacobiSNja:, y] 
is a Jacobian elliptic function. Eqn.(17) represents an 
exact solution of the Eqn.(15). The normalization con- 
dition A J Q (j>Q(z)dz = 1 can be recast from the Eqn.(17) 
as 381 



8EllipticK[M] 



EllipticK[M] - EllipticE[M] 



L 

A' 
(18) 



where EllipticE[M] = f - -^M 2 ~ ^M 3 - ...} [H 
is also another elliptic integral. Now, putting z = L/2 
and M = ^- in the Eqn.(17), we get [V^ 



L 
2 



1 EllipticK[- 
sfd VI - V 



(19) 



which can once again be recast in the form of the ground 
state energy [38| 



M = 



2ft 2 EllipticK 2 [Af](l + M) 
mL 2 



(20) 



Eqn.(19) is a transcendental equation for the maximum 
of <j>o(z). M can be determined form the transcendental 
Eqn.(18) by a graphical method 38]. Once we put the 
determined M in to the Eqn.(20), we will get the ground 
state energy. 

It is easy to check from the Eqn.(19) and (20) that, 
the allowed ranges of rj and M are < rj < ^ and 
< M < 1 respectively. In the non-interacting limit 
(M — > 0), Eqn.(20) gives the ground state energy of a 

2 j- 2 

particle in a box as /i = 2mL 2 ' ^ n ^ ne °t ner band, in 
the bulk limit (M 1), EiiipticKtM] ~ ^ ^ can ^ e P u * m 
the Eqn.(18), which together with the Eqn.(20) gives ^ 
(i = gn, which we already have used in the Eqn.(2). 

Since the equilibrium of our system is obtained by min- 
imizing the grand potential, it is relevant to get the force 
(acting on the plates), from the derivative of the grand 
potential with respect to the plate separation. For this 
process, it is necessary to keep the chemical potential un- 
changed. To get the Casimir force, one needs to subtract 
the bulk force (F m f (oo)) from the actual force (F m f(L)) 
in the Eqn.(14) with constant chemical potential, which 
can conveniently be taken as the bulk chemical potential 
fx = gn even for this finite system [40| . 

Let us now look at the Eqn.(14), where y (= 
by the definition is a constant, and a (— 2{j,m/h 2 ) is to 
be kept constant (2gnm/h 2 ) for the measurement of the 
F m f(L), and that of the Casimir force. Therefore, for 
obtaining the Casimir force, we need only to know the 77 
(in the Eqn.(14)) in terms of a and L. From the Eqn.(19), 
we can determine 77 in terms of a and L by asymptotic 
and graphical analyses. 



ASYMPTOTIC AND GRAPHICAL 
ANALYSES OF r? 

A. Asymptotic analysis 



It is easy to check from the Eqn.(19) that, for r/ — > 0, 

a — > and for 77 — > j, aL 2 — > 00. Let us define a — 

2 

a — jp. Now we can check from the Eqn.(19) that, 77 
as well as <^o(§) goes to zero as L\/E — > 0. For 77 — > 0, 
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FIG. 1: Comparison between the interpolated r\ (continuous 
line), the asymptotic 77's of the extreme limits, and the graph- 
ical solutions (dots) of r\ . The dots are obtained from the 
graphical solutions of the Eqn.(24) with the help of Wolfram 
Mathematica 5.2. The dotted, dashed and the continuous 
lines follow the Eqns.(21), (22) and (23) respectively. 



0.1 0.2 0.3 0.4 0.5' 

FIG. 2: Thick continuous line represent the right hand 
side of the Eqn.(24). The other lines represent the left 
hand side of the Eqn.(24) for x = = 1.0, 1.5, and 2. 

Intersecting points give the graphical solutions r\ = 
0.3206, 0.4264, and 0.4721 respectively. Along with other such 
points, we plot the dots of the FIG. 1. 



we can expand the right hand side of Eqn.(19) up to the 
third order in 77 and get 



n 



2 ( L\fa 



x _25 



24 V 7T 



1.04514 



Ly/a 



for < 1. 



(21) 



In the other limit, i.e. for 77 — >■ -i , Ly/a goes to in- 
finity. In this limit 77 can be expressed as 77 = i(l — 
6), where 6 — > 0. Let us now evaluate the 6 from 
the Eqn.(19), which can be expressed in the integral 

form with the substitution p = 1 — t as — v 2 L ~ = 
7=3= -J = = dp = , which says that p = 



/p( 2 -p)(l-±=£(l-2p+p 2 )) 

has a logarithmic divergence, so that most of the in- 
tegral would come from p — > 0. In this limit the above 
integral can be approximated with the first order in 5 



and p, as ^ 
asymptotically gives 



Jo 



dp 



y/MS+p) 



V^sinh" 1 ^), which 



77 = -cotlT 
' 2 



for 



Ly/a 



> 1. 



(22) 



Let us now find a smooth function for 77 as an interpo- 
lation for the whole range < i] < |, in such a way, 
that, it fits to the extreme ends of 77. Since 4>q(z) = 
\J~^t&ah.{z\Ja/2) is a solution of the Eqn.(10) for an un- 
bounded situation, we can take a trial for the bounded 

system that, <po(j) would be close to A/ftanh(- 



/a/2\ 



With this consideration, we can take a trial interpolation 
as g3 



1 



-tanh 




for 0< 



< 00. 



(23) 



For ^5 > 1 ; the trial interpolation in the Eqn.(23), fits 

obviously well with the 77 of the Eqn.(22). For < 1, 
the trial interpolation fits also well with the 77 of the 
Eqn.(21). We can see how well the interpolation fits 
with the two asymptotic 77s in the FIG. 1. Although 
the asymptotic analysis is easier to express mathemati- 
cally, yet the interpolation is not beyond the doubt for 
the region 0.5 < x < 3 in the FIG 1. For this reason, we 
need the graphical analysis to support the interpolation 
formula. 



B. Graphical analysis 



Eqn.(19) can once again be recast in terms of x — 



as 



Ixji - 7/ = EllipticK 



1-7, 



(24) 



For a given x, we plot the left hand side of the Eqn.(24) 
with respect to 77 in the FIG. 2. In the same figure, we 
also plot the right hand side of the Eqn.(24) with respect 
to 77. The intersection point is the graphical solution of 
77 for the given x. For a given set {x}, we can get the 
corresponding graphical solution {77}. We plot this set of 
points in the FIG. 1, and can compare the interpolation 
formula (Eqn.(23)) with the graphical solutions of 77. 

Now, we see in the FIG. 1 that, for the entire region 
< x < 10, the interpolation fits very well with the 
graphical solution of 77. So, for the rests of this paper, 
we will calculate the Casimir force from our interpolation 
formula (Eqn.(23)). 
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VI. CASIMIR FORCE 
A. Mean field part of the Casimir force on a BEC 

Since in the bulk limit tj — > |, we get the bulk force 

from the Eqn.(14) as F m/ (oo) = . Subtracting 

this bulk force from the mean field force in the Eqn.(14), 
we get the mean field Casimir force 

„ , , NAh 2 a 2 fl A 
C mf (L,n) = —-^--nil-n)], (25) 

which is obviously attractive, and can be simplified with 
the Eqn.(23) in the following form 



C m f{L, n) 



NAti 2 



imb 



-sech 



L 2 c 



8 



(26) 



C mf (L,p) /C qf (L,p) , <9(n c a 3 s ,p) 



2.5 




10 12 14 



FIG. 3: The continuous line follows the Eqn.(32), and repre- 
sents the total Casimir force in units of \Cp hn \{n c a i 3 )~ 1 ' 2 . The 
dotted line represents the ratio of the mean field (Eqn.(28)) 
and fluctuations (Eqn.(31)) parts of the Casimir force in units 
of 10 3 . Here we set L = 10" 7 m and n c ai = 4.0117 x 10 -5 . 



Now, putting the values of the constants a, b and g, we 
recast the Eqn.(26) as 



C m f(L, n) 



A2nh 2 a,n 2 



-sech 



m 



8nL 2 a s n — ir 2 



(27) 



It is clear from the Eqn.(27) that, for a s > 0, there exists 



a lower critical density n c 



of atoms, at and below 



" c ~ SL 2 a 3 

which we can not write p = gn. The critical density also 
defines a critical phonon speed v c = \J ' gn c /m. Now, 
the Eqn.(27) can be expressed in terms of the unit less 
'density' p = -7- as 



c mf (L, P ) = -|q 



phn I 



,3^-1/2 



x p^sech \ -z-(p~ 1) 



(28) 



where Cp hn = — ^ol 4 " is the primary quantum Casimir 
force (Eqn.(7)) for the critical density. 



B. Quantum fluctuations part of the Casimir force 
on a BEC 



Now, replacing the sum over pj_ by an integration from 
zero to an ultraviolet cutoff (A), we recast the Eqn.(29) 

as 



j(£,n) 



1 y-v 2nA 

2^(2^ 

yl 4 - (nhj/ L) 4 

A^/Amgn 
A' - (nhj/ L) 7 
&96{mgn) 2 



A 3 - (Trhj/Lf 



A 5 - (nhj/L) 5 
AOmgn 
A 2 



(30) 



The Casimir force {C q f{L,n) = — J^Eo{L,n)) from the 
low obtaini 

i^'vi i) + (^) 4 C(-4) _ ::.-/■ )\i 



Eqn.(30) is now obtained as C qf (L, n) = - 



L 4 1 ^4mgnL 5 8mgnL 6 128(mgn) 2 L s "')"> 

which can be recast (with the regularized zeta) in terms 
of p as 



C gf (L,p) = -\C c phn \p 1 ^[l--^- 



_5_ 

42p 



1 



64p 2 



(31) 



The first term of the Eqn.(31) is the dominating term 
for the quantum fluctuations part of the Casimir force as 
already obtained in the Eqn.(7) and in the Refs.[25l [35|. 



Let us now go back to the Eqn.(6) whose second term 
(eo(L,n)) is relevant for quantum vacuum fluctuations 
as well as for the quantum Casimir force. Expanding 

e(p±,j) in terms of p = \jy\ + ■ 7r2 ^ 3 " , we can write the 
second term as 



e Q (L,n) = 



V 



yJAmgn &mgn 



I28m 2 g 2 n 2 



gn 



(29) 



C. Total Casimir force on a BEC 

The total Casimir force (C(L,p)) is the sum of the 
right hand sides of the Eqns.(28) and (31), and can be 
given in units of |Cp ?ln |(n c ag) _1 / 2 by 



V 2 

+ ^^)p l ' 2 (l 



1 



42p 64p 2 



(32) 
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where n c a 3 is the parameter which determines the di- 
luteness of the condensate. Diluteness is one of the nec- 
essary conditions for achieving BEC. For the diluteness 
we must have rt c a 3 <§C 1. For 23 Na atoms (a s = 19.1ao 
[H]), and for L = 1CT 7 m, we have n c = 3.88 53 x 10 22 /m 3 
and n c a 3 = 4.0117 x 10 -5 . It is to be mentioned that, 
n c = 3.88 53 x 10 22 /m 3 is achievable in the 3D harmonic 
traps of the Bose-Einstein condensation experiments [43[ . 
For n c a 3 = 4.0117 x 10~ 5 , we plot the right hand side of 
the Eqn.(32) in the FIG. 3 with respect to the 'unitless 
density'. 

From the Eqn.(32), we get the total Casimir force 
as C(L,p) = \C c phn \{n c al)^ 2 x d{n c al,p). With the 
above parameters, and with A — 10~ 6 m 2 , we have 
\C^ hn \ = 1.318626 x 10~ 15 N, and C(10- 7 m, 1) = 
-2.7687 x 10~ 12 N which is certainly a measurable 0,11] 
quantity. 

From the Eqns.(28) and (31) we can also compare the 
mean field and quantum fluctuations' contributions to 
the Casimir force. We also plot their ratio (in units 
of 10 3 ) in the FIG. 3. For the above parameters, the 
mean field part is 2425 times stronger than quantum 
fluctuations part at p = 1. However, at higher density 
(p > 10.96) the fluctuation part dominates over the mean 
field part. 

VII. CONCLUSIONS 

This paper describes the theoretical calculations on 
the Casimir force (as well as pressure) exerted by a self- 
interacting Bose-Einstein condensate on planar bound- 
aries that confine the condensate. The calculations have 
been performed within the quadratic fluctuations over 
mean field level. 

The mean field part of the Casimir force is the conse- 
quence of the inhomogeneity of the condensate. For high 
density p > 10, the condensate between the two plates 
becomes essentially homogeneous, and the fluctuations 
part dominates over the mean field part. 

All our calculations or results are valid only for the re- 
pulsive interaction (a s > 0). For a s < 0, the condensate 
becomes unstable beyond a critical number of particles 
[lif . For T — > 0, Casimir force (in the Eqns.(7) and 
(27)) becomes zero [H, [2j| only in the non-interacting 
limit (a s — > 0), otherwise it is nonzero, attractive and 
experimentally measurable. 

For the calculation of the Casimir force, it is necessary 
to change the system size at constant chemical poten- 
tial by allowing particle exchange between the conden- 
sate and a particle reservoir. This is certainly convenient 
for the calculation, but difficult for an experiment. The 
Casimir force can be measured in this way only above 
the critical density of the atoms. 

The nature of the Casimir force on the BEC (as shown 



in the FIG. 3), is similar to that observed [l5| for the 4 He 
film above its lambda point. 

For the consideration that the condensate is confined 
between two infinitely large parallel plates along the x—y 
plane, the equation of motion of the condensate is rele- 
vant only in the z direction. For the same reason, we 
have solved the ID nonlinear equation (Eqn.(8)). But, 
practically the plates are of finite area. If A be small, the 
equation of motion of the condensate would be relevant 
in the x, y and z directions. Calculation of the Casimir 
force for the plates of small area would of course be an 
interesting problem. In this situation, we need to solve 
the Eqn.(8) by replacing <fi(z) by 4>(r). This replacement 
would make the Eqn.(8) a 3D nonlinear equation, which 
is very difficult to solve due to the fact that separation 
of variable technique can not be applied to the nonlin- 
ear problems. Hence, the 3D calculation for the same 
problem is difficult but interesting. 

For a realistic case, 3D harmonically trapped BEC is 
to be placed between two plates. The condensate would 
exert a force on the plates if the plates are kept a little 
away from the condensate. For this realistic geometry, 
the magnetic field may shield the condensate from being 
destroyed by the (Van der Waals) interaction between the 
plates and the particles. If the plates are kept very far 
away from the condensate, then the (bulk) force acting 
on the plates would be zero. So, for the realistic geom- 
etry, the force acting on the plates itself is the Casimir 
force. The experimentalists [20| already investigated the 
Casimir-Polder effect by putting a single plate away from 
a condensate. Hence, we may expect that the experimen- 
talists might be able to measure the Casimir force for a 
realistic geometry of BEC. 

However, we could not get the Casimir force on a 3D 
harmonically trapped interacting BEC. How to get this is 
an open question. How to calculate the Casimir force on 
an inhomogeneous interacting Bose gas for T > 0, is also 
an open question. The same problem with interacting 
fermions may also be an interesting issue. 
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